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摘要 :大 功率 中 高 速 柴油 机 是 众多 关键 装备 的 核心 动力 设备 ， 在 舰 船 、 核 电 等 重要 行 
业 领 域 有 着 广阔 的 应 用 。 由 于 复杂 的 机 械 结构 和 多 变 的 工作 环境 ， 柴 油 机 设备 事故 时 有 发 
生 ， 机 械 故 障 和 设备 损毁 等 问题 也 日 益 凸 显 ， 故 障 频 发 成 为 这 些 设备 运行 的 一 大 难题 。 柴 
油 机 恶性 故障 的 发 生 经 常会 导致 严重 的 经 济 损失 和 人 员 伤 亡 ， 因 此 ， 采 用 合理 的 方法 获取 
故障 特征 ， 研 究 对 柴油 发 动机 的 故障 检测 和 早期 诊断 方法 ， 实 现 对 柴油 发 动机 的 故障 预警 
并 及 时 诊断 ， 对 提高 大 功率 中 高 速 柴油 机 运行 的 稳定 性 、 设 备 的 安全 性 和 企业 的 经 济 性 有 
着 重要 的 作用 与 价值 。 

本 文 提 取 柴 油 机 角 域 、 时 域 、 频 域 22 种 振动 特征 组 成 特征 集 ， 针 对 柴油 机 应 用 场合 众 
多 且 工 况 复 杂 的 特点 ， 在 仅 有 少量 工 况 样 本 的 情况 下 ， 本 文 基于 神经 网 络 实现 工 况 与 马 氏 
空间 关键 参数 的 映射 关系 ， 构 建 变 工 况 下 柴油 发 动机 的 健康 基线 。 
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Abstract: Many key equipment rely heavily on the high-powered medium to high 
speed diesel engine, which has a broad range of applications in areas such as 
ships and nuclear power. Due to the complex mechanical structure and ever- 


changing working environment, diesel engine equipment accidents occur frequently, 


and mechanical failures and equipment da 
Frequent failures have become a major 
equipment. The occurrence of malignant 
serious economic losses and casualties. 
obtain fault characteristics, researchi 
methods for diesel engines, achieving 


diesel engine faults, plays an important 


of high-power, medium and high-speed diesel engine operation, 


and enterprise economy. 


age are becoming increasingly prominent. 


challenge for the operation of these 
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Therefore, using reasonable methods to 
ng fault detection and early diagnosis 


early warning and timely diagnosis of 
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equipment safety, 


This article extracts 22 vibration characteristics of diesel engines in the 


angular domain, time domain, 


and frequency domain to form a feature set. 


In 


response to the diverse and complex operating conditions of diesel engines, with 


only a small number of operating samples, 


map the relationship between operating conditions and key parameters 


this article uses neural networks to 
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Markov space, and constructs a healthy baseline for diesel engines under variable 
operating conditions. 
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1 绪论 


国外 对 于 复杂 工 况 条 件 下 的 故障 诊断 起 步 较 早 。 早 在 1999 年 ， Raj K 就 发 现 
采用 模糊 ARTmap《〈 自 适应 共振 理论 ) 神经 网 络 非 常 适合 解决 各 种 系统 和 故障 条 
件 下 的 复杂 故障 分 类 问题 。 他 重点 研究 了 应 用 于 特定 问题 的 人 工 智能 技术 ， 然 后 
描述 用 于 训练 Fuzzy ARTmap 神经 网 络 的 方法 ， 此 外 ， 将 该 分 类 技术 与 基于 误差 
反 向 传播 (EBP) 训练 算法 的 神经 网 络 (NN) 技术 进行 了 比较 ,结果 表明 ， 前 者 
更 适合 解决 复杂 多 回路 传动 系统 的 故障 诊断 问题 息 。2001 年 , 日 本 的 Toshio 
Toyota 针对 提取 最 重要 的 兆 参 数 的 函数 困难 的 问题 ， 提 出 了 一 种 基于 遗传 算法 
的 “征兆 参数 自动 函数 生成 ”方法 , 并 将 该 方法 应 用 于 真实 的 机 械 故障 诊断 问题 ， 
该 方法 能 快速 生成 关键 征兆 参数 函数 申 。 在 2007 年 第 52 届 ASME Turbo Expo 
上 ，Igor Loboda 针对 燃气 涡轮 机 设备 和 发 动机 的 运行 条 件 〈 控 制 变量 和 环境 条 
件 ) 变化 很 大 ， 提 出 利用 热力 学 模型 来 模拟 气 路 故障 ， 并 使 用 神经 网 络 来 进行 故 
章 定 位 ， 重 复 执 行 该 技术 并 记录 诊断 ， 基 于 这 些 诊 断 和 预先 已 知 的 故障 ， 然 后 计 
算 正 确 的 诊断 概率 , 分 析 了 运行 条 件 对 诊断 过 程 的 影响 , 提出 模拟 在 可 变 操作 条 
件 下 获得 的 测量 值 的 诊断 处 理 的 方法 辐 。M. Ben Khader Bouzid 提出 了 一 种 基于 
使 用 前 馈 多 层 感知 人 工 神经 网 络 (ANN ) 的 永 磁 同步 发 电机 定子 故障 自动 智能 诊 
断 系 统 , 通过 反 回 传播 训练 算法 执行 ,并 通过 对 不 同 定子 故障 和 不 同 运行 工 况 下 
电机 负 序 电压 行为 的 深入 研究 和 分 析 , 选择 了 两 个 新 的 故障 鲁 棱 性 指标 作为 神经 
网 络 的 输入 四 。Shahid 等 人 使 用 磁性 感应 器 来 监控 转速 ， 并 将 所 测 得 的 信号 转化 
为 角 域 信号 , 然后 将 负载 按 比 例 划分 为 5 个 等 级 ,利用 多 层次 感知 机 实现 了 对 发 
动机 负载 的 实施 分 类 四。 

相 比 之 下 , 尽管 中 国 在 此 领域 的 发 展 落后 于 其 他 发 达 国 家 , 但 已 经 取得 的 较 
为 突出 的 进展 。 柴 油 机 的 工作 状态 与 燃烧 产生 的 功 有 一 定 的 关联 ， 因 此 ， 通 过 对 
其 工作 条 件 的 分 析 ， 可 为 其 优化 设计 提供 依据 ， 从 而 提高 柴油 机 的 工作 效率 。 国 
内 自行 研发 的 故障 监测 系统 发 展 迅 速 , 已 经 积累 了 大 量 宝贵 的 数据 和 经 验 , 并 且 
广泛 应 用 于 电力 、 石油 和 化 工 等 领域 。 但是, 在 这 方面 的 技术 研发 和 应 用 还 远 远 
落后 于 发 达 国 家 ， 需 要 继续 加 大 研发 力度 ， 不断 创 新 与 实践 。 目 前 ,国内 许多 学 
者 都 在 研究 中 采用 基于 煤油 机 和 氏 压 辨识 的 方法 来 判断 煤油 机 的 工作 状态 。 李 柱 国 


和 其 他 研究 人 员 提 出 了 一 种 新 的 故障 预警 方法 中, 该 方法 基于 包 盖 振动 信号 包 络 
线 。 刘 建 敏 等 人 利用 缸 盖 在 稳定 状态 下 的 振动 信号 ,利用 遗传 算法 和 神经 网 络 相 
结合 的 方法 实现 柴油 发 动机 缸 内 压力 识别 与 判断 004。 常 春 等 人 基于 时 域 - 频 域 相 
关 分 析 法 ， 从 发 动机 的 燃烧 室 中 提取 出 发 动机 的 燃烧 性 能 指标 ， 然 后 基于 径 向 基 
神经 网 络 来 估算 发 动机 的 内 压 00。 相 对 于 气缸 压力 而 言 ,， 柴油 机 的 转速 和 负载 是 
更 能 直接 反映 柴油 机 工作 状态 的 指标 。 通 过 对 上 荣 油 发 动机 振动 信号 的 分 析 ， 可 以 
直接 提供 柴油 机 控制 所 需 的 工 况 信息 , 还 可 以 了 解 某 油 机 工作 状态 对 振动 信号 的 
影响 。 同 时 ， 这 种 方法 也 可 以 应 用 于 上 某 油 机 的 振动 信号 的 识别 与 诊断 等 方面 。 江 
志 农 等 人 利用 多 种 信号 特性 ， 包 括 温 度 、 压 力 、 振 动 等 ， 构 建 了 一 种 新 的 t- 分 布 
区 域 艇 入 方法 (t-distributed stochastic neighbor embedding) ， 实 现 对 数据 的 有 效 
降 维 , 并 且 进 一 步 研究 了 一 种 新 的 k 最 近邻 分 类 算法 ,实现 了 对 柴油 机 运行 情况 
的 准确 分 类 02。 

传统 的 健康 基线 构建 方法 通常 忽略 了 多 变 因素 的 影响 , 而 依靠 经 验 或 者 故障 
机 理 。 对 于 柴油 机 这 种 结构 复杂 ， 故 障 种 类 较 多 的 设备 来 说 ， 由 于 缺乏 判断 经 验 
和 数据 ,使 用 传统 方法 来 构建 健康 基线 往往 会 面临 较 大 的 困难 。 华 成 等 人 提出 了 
一 种 新 的 实时 可 靠 性 评估 方法 , 该 方法 基于 故障 诊断 和 监测 技术 , 通过 Parzen 窗 
非 参 数 概率 密度 估计 法 ,构建 滑动 概率 神经 网 络 ， 不断 滑动 时 间 窗 ， 从 而 实现 对 
柴油 机 的 实时 可 靠 性 评估 M1, 这 一 方法 比 传统 的 经 验方 法 更 加 准确 , 而 且 能 够 有 
效 地 解决 实际 情况 与 假设 不 符 的 问题 。 殖 玉 良 等 人 利用 滑动 概率 神经 网 络 , 建立 
退化 基线 ， 并 基于 ARMA 模型 对 关键 参数 进行 预测 中。 大 洛 帆 利用 BP 神经 网 
络 模型 构建 出 具有 较 高 容错 性 和 良好 泛 化 性 的 柴油 机 故障 诊断 系统 , 以 有 效 地 检 
测 和 监控 某 油 机 零 部 件 的 磨损 情况 051。 

可 以 看 出 ,对 于 柴油 机 故障 的 发 生 进行 预测 和 评估 ,实现 健康 基线 的 构建 是 
一 种 十 分 有 效 的 方法 ; 而 针对 柴油 机 系统 复杂 以 及 工 况 多 变 的 问题 ， 基 于 神经 网 
络 对 柴油 机 实时 监测 和 故障 预测 有 着 很 好 的 作用 与 效果 09。 

本 文 的 研究 对 象 是 大 功率 中 高 速 柴油 发 动机 , 研究 目标 是 构建 变 工 况 下 大 功 
率 中 高 速 柴 油 发 动机 的 健康 基线 ， 实 现 发 动机 的 故障 预警 。 本 文 引入 11 种 不 同 
工 况 下 煤油 发 动机 的 振动 信号 ,提取 信号 的 特征 值 ， 构 建 马 氏 距 离 空间 , 针对 复 
杂 多 变 的 工 况 ,基于 深度 学 习 神经 网 络 来 对 变 工 况 下 的 数据 进行 拟 合 , 对 数据 进 
行 训练 ， 最 终 构 建 变 工 况 下 的 健康 基线 。 
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2 ” 变 工 况 下 振动 信号 特征 数据 集 的 分 析 与 处 理 


2.1 变 工 况 下 柴油 机 振动 信号 特征 的 分 析 
本 文 以 11 种 大 功率 中 高 速 工 况 为 例 ， 其 转速 和 负载 如 表 2-1 所 示 。 
表 1 11 种 工 况 对 应 转速 、 负 载 表 


工 况 数 1 2 3 4 5 6 7 8 9 10 11 
转速 / 900 950 1000 1000 1000 1050 1050 1050 1100 1100 1100 
rpm 
负载 / 0 0 0 100 200 0 100 200 0 100 200 
N*m 


分 别提 取 每 组 样本 五 段 冲击 区 间 的 开启 相位 、 关 闭 相位 和 区 间 宽 度 3 种 特 
征 、15 种 时 域 特征 、4 种 频 域 特征 ， 并 按 顺 序 排列 起 来 ， 组 成 具有 5*22 个 特征 
的 数据 集 矩 阵 。 按 同样 的 方法 求 出 11 种 工 况 下 300 组 样本 信号 , 得 到 11*300*S*22 
的 信号 特征 集 矩 阵 。 

将 11 种 工 况 下 五 段 冲击 的 相位 区 间 整 理 为 图 2-3 如 下 : 
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图 1 11 种 工 况 下 的 五 段 冲 击 的 相位 开 闭 区 间 
其 中 横 坐 标 从 左 到 右 依次 为 $ 个 区 间 开 始 、 关 闭 相 位 ,其 两 个 相 邻 相位 的 距 
离 即 为 区 间 宽 度 ， 纵 坐标 从 下 往 上 依次 为 11 种 工 况 (用 红线 隔 开 )。 
分 别 求 出 每 组 的 22 种 特征 ， 组 合成 11*300*5*22 种 特征 ， 得 到 如 下 图 2-4。 
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图 2 11 种 工 况 下 300 组 数据 的 22 个 特征 集 

其 中 x 轴 0:110 表示 5 个 区 间 的 22 个 特征 ， 共 110 个 特征 ，y 轴 0:3300 表 
示 11 种 工 况 下 每 种 工 况 的 300 组 数据 ， 共 3300 组 数据 。 
2.2 变 工 况 下 柴油 机 振动 信号 特征 的 处 理 

在 提取 出 柴油 机 正常 样本 的 特征 集 后 ， 需 要 对 数据 进行 规范 化 度量 ， 以 便于 
更 好 的 对 柴油 机 系统 进行 健康 评估 。 由 于 各 种 特征 的 量 纲 不 尽 相 同 ， 且 样本 与 样 
本 之 间 有 着 一 定 的 联系 ,， 故 本 文选 择 建 立马 氏 参 考 空间 , 通过 求 出 样本 点 与 样本 
特征 之 间 的 马 氏 距离 的 方式 ， 对 样本 数据 进行 处 理 和 来 实现 健康 基线 的 构建 。 

本 文 马 氏 参考 空间 的 数据 是 由 正常 状态 下 柴油 发 动机 样本 的 特征 来 确定 的 ， 
即 第 2.1 节 求 出 柴油 机 300 组 样本 信号 5 段 冲 击 区 间 的 22 个 特征 所 组 成 的 110 
个 特征 ， 得 到 的 300*5*22 的 信号 特征 集 矩 阵 。 

在 构建 马 氏 空间 函数 后 ， 代 入 数据 ， 分 别 求 出 构建 出 的 马 氏 空间 的 均值 ， 方 
差 和 协 方差 矩阵 ， 最 后 求 出 正常 样本 的 300 组 马 氏 距离 (MD) 如 图 所 示 。 
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图 3 单 工 况 下 的 马 氏 距离 (MD) 
3. 变 工 况 下 柴油 发 动机 健康 基线 的 构建 


柴油 发 动机 在 各 个 领域 都 得 到 了 广泛 的 应 用 , 但 是 它们 在 不 同 的 环境 下 的 运 
行 转速 和 负载 要 求 也 各 不 相同 。 这 时 ， 需 要 构建 一 种 变 工 况 的 柴油 发 动机 的 健康 
基线 ， 以 适应 多 变 的 场合 。 

以 工 况 信息 [转速 、 负 载 ] 的 二 维 同 量 为 输入 ， 如 [1100，200] 代 表 转 速 为 
1100rpm， 负 载 为 200Nm, 分 别 以 均值 向 量 X、 标 准 差 s 和 相关 算 阵 C 为 输出 ， 训 
练 解码 器 网 络 ， 完 成 映射 关系 的 拟 合 。 将 生成 的 均值 向 量 均值 向 量 X、 标 准 差 s 和 
相关 算 阵 C 代入 到 马 氏 距离 相关 公式 中 ， 即 可 求 出 该 工 况 下 的 马 氏 距离 。 

3.1 神经 网 络 的 构建 

本 文 利用 Python 的 Keras 库 来 构建 深度 学 习 模 型 ， 通 过 神经 网 络 对 数据 进 
行 拟 合 ， 其 步 又 如 下 : 

(1) 构建 一 个 包含 输入 层 、 编 码 层 和 解码 层 的 自 编 码 器 模型 ， 用 于 将 输入 
数据 进行 降 维 和 重 构 。 其 中 编码 层 和 人 解码 层 都 包含 一 个 隐藏 层 , 通过 训练 自动 调 
整 权重 参数 ， 从 而 实现 对 输入 数据 的 有 效 编码 和 解码 。 

通过 采用 自 编码 器 ， 将 输入 数据 进行 降 维 和 重 构 ， 将 其 组 织 成 55 行 3 列 的 
矩阵， 每 组 数据 包含 $ 个 冲击 区 间 ， 总 计 11 组 ， 每 组 包含 55 个 样本 ， 第 一 列 代 
表 原 始 数值 ， 第 二 列 和 第 三 列 分 别 指示 11 种 不 同 工 况 下 的 转速 和 负载 。 


(2) 将 第 2 章 中 的 11 种 工 况 下 的 数据 进行 代入 ， 并 将 11*110 的 特征 集 逢 
阵 重 新 排列 ， 将 其 中 的 每 一 行 中 的 22 个 特征 的 冲击 区 间 分 为 5 组 ， 组 成 55*22 
的 矩阵 ， 为 了 方便 数据 的 处 理 ， 对 其 进行 归 一 化 处 理 。 同 时 ,记录 每 一 列 的 最 大 
值 和 最 小 值 ， 以 便 后 续 进 行 反 归 一 化 操作 。 

(3) 建立 一 个 有 5 层 的 解码 器 的 模型 ， 具体 为 一 个 有 着 三 个 特征 的 输入 层 
定义 了 一 个 30 个 神经 元 的 全 连接 层 , 一 个 100 个 神经 元 的 全 连接 层 , 一 个 50 个 
神经 元 的 全 连接 层 和 一 个 22 个 神经 元 的 全 连接 层 。 这 样 ， 就 构建 了 一 个 有 3 个 
输入 节点 、22 个 输出 节点 的 解码 器 模型 。 

(4) 设 定 模型 的 训练 参数 迭代 次 数 epochs 为 1000 次 ， 学 习 率 Ir=0.001。 

(5) 通过 decoder 解码 器 来 训练 并 预测 模型 ， 并 利用 Adam 优化 器 以 及 均 方 
误差 (MSE) 来 优化 其 表现 ， 从 而 实现 更 高 效 的 数据 人 处理 。 

神经 网 络 结构 示意 图 如 下 图 所 示 : 
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图 4 神经 网 络 解码 器 示意 图 


3. 2 拟 合 效果 

下 图 为 11 种 工 况 5 个 冲击 的 22 个 特征 的 X 向 量 〈 均 值 向 量 ) 拟 合 效果 ， 左 
图 为 实际 特征 ， 右 侧 为 拟 合 出 来 的 特征 ,将 上 述 特征 按 列 归 一 化 后 的 拟 合 效果 如 
下 图 。 


归 一 化 均值 向 量 
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图 $ 神经 网 络 训练 前 后 的 (均值 向 量 ) 对 比 
通过 拟 合 前 (图 5 左 ) 与 拟 合 后 (图 5 右 ) 对 比 ， 并 将 原始 输入 与 重 构 输出 
进行 比较 ， 以 评估 模型 重 构 输入 的 准确 性 ， 求 得 均 方 误差 MSE) 为 0.000360， 
可 以 看 出 拟 合 效 果 相 对 较 好 。 
相应 的 将 s 向 量 ( 标 准 差 向 量 ) 和 相关 和 矩阵 C 进行 拟 合 。 


归 一 化 均值 向 量 
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图 6 神经 网 络 训练 前 后 的 s( 标 准 差 ) 对 比 
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图 7 神经 网 络 训练 前 后 的 C( 相 关 和 矩阵 ) 对 比 


通过 绘制 训练 损失 随时 间 变 化 的 图 表 ， 也 可 以 显示 损失 随时 间 不 断 减 小 。 
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图 8 神经 网 络 训练 前 后 的 C( 相 关 和 矩阵 ) 对 比 


将 拟 合 后 的 数据 代入 马 氏 距离 公式 ， 构 建 11 种 工 况 下 的 马 氏 空间 ， 求 得 拟 


合 后 的 马 氏 距离 如 图 4-6 所 示 。 
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图 9 训练 后 数据 的 马 氏 距离 
所 以 该 解码 器 网 络 完成 了 工 况 信 息 到 [均值 向 量 x、 标 准 差 s 和 相关 和 矩阵 C] 的 
拟 合 ， 即 只 需要 输入 任意 工 况 信息 [转速 ， 负 载 ] 二 维 向 量 ， 即 可 输出 一 组 对 应 的 
[均值 向 量 x、 标 准 差 s 和 相关 和 矩阵 C]， 该 输出 用 于 构造 马 氏 空间 ， 即 可 完成 不 同 
工 况 之 间 的 对 齐 ， 实 现 变 工 况 下 的 马 氏 空间 的 构建 。 
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